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We study the effect of decoherence on a weak value measurement in a paradigm system consisting 
of a double quantum dot continuously measured by a quantum point contact. Fluctuations of the 
parameters controlling the dot state induce decoherence. We find that, for measurements longer 
than the decoherence time, weak values are always reduced within the range of the eigenvalues 
of the measured observable. For measurements at shorter timescales, the measured weak value 
strongly depends on the interplay between the decoherence dynamics of the system and the detector 
backaction. In particular, depending on the strength of the decoherence and the post-selected state, 
a stronger detector backaction might lead to an enhancement of weak values. 
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I. INTRODUCTION 

In quantum mechanics the measurement process is 
most simply described as a probabilistic event through 
the projection postulate^. While it satisfactory describes 
several simple experimental configurations, some mea- 
surement protocols, including conditional quantum mea- 
surements can lead to results that cannot be interpreted 
in terms of classical probabilities, due to the quantum 
correlations between measurements. A striking evidence 
of that is provided by the so called weak values (WV) 
obtained from the measurement scheme originally devel- 
oped by Aharonov, Albert and VaidmarP. The weak 
value measurement protocol consists of (i) initializing 
the system in a certain state | $ ) (pre-selection) , (ii) 
weakly measuring an observable A of the system by cou- 
pling it to a detector, and (iii) retaining the detector 
output only if the system is eventually measured to be 
in a chosen final state | $ ) (post-selection) . The aver- 
age signal monitored by the detector will then be pro- 
portional to the real part of the so called weak value 

$ (i)£ cak = ( $111* )/<$|*). 

The most surprising pro perty of WVs is that they can 
be complex or negative^ whereas a strong conventional 
measurement would lead to positively definite values. Af- 
ter the original debate on the meaning and significance 
of weak values^^, they proven to be a successful con- 
cept in addressing fundamental problems and paradoxes 
of quantum mechanics-^!, in defining measurements of 
controversial quantities (e.g. the time a particle spends 
under a potential barrier in a tunneling process^!) and 
counter-intuitive situations (e.g. the simult aneo us mea- 
surement of two non-commuting observable d 7 * 10 !), as well 
as in generalizing the definition of measurement!^. By 
now a number of experiments in quantum optics has re- 
ported the experimental observation of weak values and 
its application to quantum paradoxes^^. Recently a 
series of interesting works has explored the potential of 
weak values measurement protocols for precision mea- 
surements. Weak value based measurement techniques 
have been successfully employed in quantum optics ex- 
perimentes to access tiny effects 1 ^ and detect ulta-small 



(sub-nanometric) displacements^ 4 ^. A parallel research 
has introduced the idea of weak values also in the con- 
text of solid state state systems^HH]. Here, further works 
have shown that weak values are related to the violation 
of classical inequalities in current correlation measure- 
ments^, and a weak value measurement technique for 
ultra-sensitive charge detection has also been proposed^. 

Due to the fact that weak values stem out of quantum- 
mechanical correlations between two measurements they 
are expected to be particularly sensitive to decoherence. 
The effect of decoherence is important for weak value 
ultra-sensitive measurements where decoherence could 
suppress the amplification effect and become crucial in 
possible solid state implementations, where it is known 
that decoherence plays a significant role in most systems. 
This is in fact the case for all the actual proposed imple- 
mentations of WVs in solid state systems^ " 18 * 20 * 21 ! So 
far, the effect of decoherence on weak values has recently 
been addressed at a formal level showing how weak val- 
ues are defined in a general open quantum system 3 -!, while 
a quantitative evaluation of the effects of decoherence in 
a sp ecific system exists only for spin qubits in a simple 
Therefore, a general characterization of the 
effects of decoherence within a weak value measurement 
in an open quantum system is, besides its theoretical sig- 
nificance, a relevant step in the direction of weak value 
implementation in condensed matter systems. 

In this work we precisely address this question. We 
approach the problem by considering the effect of deco- 
herence in a paradigm system, namely a quantum point 
contact (QPC) sensing the charge in a nearby double 
quantum dot! 23 * 24 ! The model captures all the essential 
features of a continuous quantum measurement, corre- 
sponding to the typical measurement schemes of quan- 
tum states in nanoscalc solid-state systems (which is the 
case for all the above mentioned proposals), and allows 
us to fully describe the interplay between the detector 
backaction and the decoherence process. 

The key features of our analysis and the main results 
are as follows. We describe the double quantum dot as a 
two level system, \L), \R), corresponding to the electron 
being in the left or right dot. In the system dynamics, 



we introduce fluctuations of the system's parameter, e.g. 
the gate voltages, that suppress the quantum mechanical 
oscillations between these two states, at the decoherence 
time, I/7. The QPC detector, while distinguishing be- 
tween the two system states, effects the system dynam- 
ics on the scale of the backaction time, tq. As long as 
the measurements duration, tm is shorter than td, the 
measurement is weak and can lead to weak values upon 
proper post-selection. The post-selection, obtained, e.g. 
by a second detector, is effectively described as a projec- 
tive measurement on a specific state, \ipf). 

The quantity of interest in such a weak value proto- 
col is the detector signal conditional to a positive post- 
selection, with a particular attention to the apparence of 
peculiar weak values, i.e. weak values which lie beyond 
the range of eigenvalues of the measured observable. By 
taking advantage of the Bayesian formalism, which allows 
us to consider the correlations for single shot measure- 
ments, we obtain a general expression for the weak value 
in terms of the system state only — c.f. Eqs. ( 15|16 ) . 

We identify two different regimes depending on 
whether the detector readout (i.e. relaxation processes) 
are slow or fast compared to the duration of the mea- 
surement. In all cases, measurements longer than the 
decoherence time lead to the weak values bounded by the 
eigenvalue spectrum. In the former case, dubbed coher- 
ent detection, we show that the weak value is exclusively 
determined by the system dynamics undergoing decoher- 
ence — c.f. Eqs. (25). In the latter case, named contin- 
uous readout, instead, even at time scales Tm *C 175, the 
weak values are affected by t he dete ctor and the deco- 
herence dynamics — c.f. Eqs. ( 32|33 ). The interplay be- 
tween the two can enhance the corresponding weak value 
obtained for coherent detection. In particular, depending 
on the various parameters, e.g. orientation of the post- 
selection, it could give rise to weak values beyond the 
range of eigenvalues, whereas a coherent detection would 
not — c.f. Fig. [5] 

The paper is primarily separated into four parts. In 
Section [n] we present the model and its description in 
terms of the Bayesian formalism. Hereafter, section |III| 
presents a general expression for the weak value in terms 
of the density matrix of the combined qubit-detector sys- 



(a) 



tern. In section IV we discuss the system dynamics of the 
two regimes of "coherent detection" (IV A) and contin- 



uous readout (IV B). In the former regime, the detector 



and system dynamics are decoupled from each other in 
the weak measurement regime, in the latter regime, we 
show that the detector induced decoherence and the in- 
trinsic system decoherence act together in affecting the 
measurement outcome. Section [V] summarizes our re- 
sults. 



II. THE MODEL 

The system under study is illustrated in Fig. [T] It 
consists of a double quantum dot with an adjacent QPC, 
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FIG. 1. (a) Scheme of a double quantum dot with a nearby 
QPC measuring its charge states, labelled as \L), \R). (b,c) 
Scheme of the detection mechanism. The conductance of the 
detector is sensitive to the qubit's occupancy: the tunneling 
amplitude in the QPC is fi+iJfi when the right dot is occupied 
(b) and fl — SQ, when the left dot is occupied (c). 



which serves as a detector continuously measuring the 
charge state in the double dot. 

We consider the case where, due to charging effects, 
the double dot can host only one electron in two orbital 
levels, \L) and \K) corresponding to the lowest orbital 
level for the left and right dot respectively. In this case 
the double dot can be thought as a charge qubit. In the 
entire paper we will refer to the "dot plus QPC" as the 
combined qubit-detector system which is then described 
by the Hamiltonian 



H(t) = H (t) + H lnt + H, 



(1) 



where H (t) denotes the Hamiltonian of the double dot, 
-ffint characterizes the system-detector interaction and H 
describes the Hamiltonian of the detector. Specifically 



H (t) = e(t)a z + A(t) a x 

H + H int = -( AH+ + AH~ ® a. 



(2) 
(3) 



where 



AH + = ^2 E r a\a r + Eiajai + ^ ( Q.i r a\a r + h 



AH =2 ^ \ SQi r ala r + h.c. 



(4) 
(5) 



Here a z = |R)(R| - |L)(L|, <j m = | R>< L | + | L >< R | 
and I R) , | L) are, without loss of generality, the eigen- 
states of the operator measured by the QPCpS The on- 
site energy difference, 2e(i), and the tunneling between 
such states, A(t), are controlled by the externally applied 
voltage biases, Pl,Pr, Vk- In Eqs. (3j5l, the operators 
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in the QPC-space are indicated by a and we will do 
so henceforth. The creation and annihilation operator 
in the two leads of the QPC are denoted by a] and Oj 
respectively, where i = l,r refer to the left and the right 
reservoir, respectively. B, characterizes the energy of the 
reservoir states, which are maintained at the correspond- 
ing Fermi-energies pl — PR + eVd > Pr, and fij )t . ± 5Cli <r 
denotes the tunneling amplitude between the reservoirs 
when | R) or | L) is occupied. 

This mod el has been extensively studied in the lit— 
erature^^^H. The current through the detector di- 
rectly measures the "position" of the electron in the 
dot. Under the assumptions of uniform tunneling ma- 
trix elements, SQi_ r = 50,, Q,i jT = ft, and density of 
states in the QPC's reservoirs, pz, Pr, the average cur- 
rent for an electron being in the left or right dot, reads 
(7)± = eD/2 ■ (1 ± 5 ft /ft) 2 , where D = TV d /{2n) = 
2irpip r n 2 Vd- Throughout the work we set H = 1. The 
QPC's shot noise power Si — el(l — T) sets the timescale 
t d ~ S I /{I+-I-) 2 « l/D(ft/5ft) 2 (l-T) needed to dis- 
tinguish the detector's signal form the background noise. 
The weak measurement regime is then identified by mea- 
surements of duration tm <C td, which can be con- 
trolled in principle by tuning the system-detector cou- 
pling, the duration of the measurement, or the voltage 
bias across the QPC. It is important to notice that, in 
the present model of the QPC, we assume that the equi- 
librium state of the reservoirs is not affected by fli r . In 
other words the backaction of the detector is only due to 
out-of-equilibrium effects (V^ = 0), while we are neglect- 
ing possible equilibrium backactions, e.g. orthogonality 
catastrophe dephasing! 29 ! 30 ! 

An important aspect of the model is that the Hamil- 
tonian in Eq. ([IJ describes the qubit and the detector as 
a closed quantum system. However, the QPC is contin- 
uously converting the information about the state of the 
system in a classical — macroscopic — information out- 
put which is the current. In other words, while Eq. ([!]) 
will evolve the detector to a coherent superposition of 
states with different charges in the reservoirs, the classi- 
cal knowledge of the current would correspond to a well 
defined number of electrons in each reservoir. A solution 
of this problem, as pointed out in Ref. [28j consists in 
introducing a macroscopic pointer which interacts with 
the detector. The pointer provides in fact an effective 
description of the various relaxation processes confining 
the QPC electrons in one of the two QPC reservoirs. The 
pointer can be modelled to interact instantaneously at 
certain times, t%, t% . . . t& . . . tjy with to = and tjy = r, 
and reads out the change of the number of electrons in the 
right reservoir, pR, as schematically depicted in Fig. [2] 
At any time t — tk the pointer reads the number of elec- 
trons, rrik, transmitted to the collector within the time 
interval Atk = tk+i — tk, and collapses the qubit-detector 
system onto a corresponding state depending on the mea- 
sured value rrik . Note that the introduction of the pointer 
results in a new timescale Atk in the problem which is a 
free parameter in our model. We will discuss in the re- 
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FIG. 2. (a) Qubit plus detector interaction with a classical 
pointer reading out the number rrik of electrons having pen- 
etrated to the detector's right reservoir at certain discrete 
times t — tk- Illustration of the continuous readout (b) and 
coherent (c) measurement regimes, respectively. 



suits the different regimes corresponding to the relation 
of this timescale to other timescales in the problem. 

The model allows us to include decoherence by consid- 
ering fluctuations, which naturally arise in the system, 
of the external parameters, namely the voltage biases 
Vl,Vr, V/j. In turn they lead to fluctuations of the dots' 
parameters, e(t) = e + and A(t) — A + n(t) 

around their average values eo and Ao- Since the ex- 
plicit relation between V^, Vl, Vr and eo, A depends 
on microscopic details, we effectively assume here a gen- 
eral linear relation and the effect of decoherence is there- 
fore described by replacing the dot's Hamiltonian by 
H (t) H + He(t), with 



H = e a z + A a x = wm • a, 



(6) 
(7) 



Here u — e 2 + A 2 , is the oscillation frequency of the 
system and m defines the corresponding eigenstates. To 
this regard, the index i labels the different independent 
decoherence sources, ideally corresponding to the inde- 
pendent voltage sources. For each of them k t indicates 
the direction of the fluctuations with |kj| = 1 and 
is assumed to be a Gaussian white noise, i.e. 

(&(s)) e = 0, ((i(s 1 )^ j (s 2 )) $ = li S(s 1 -s 2 )S ij (8) 

where describes the strength of the correlation func- 
tion. For the sake of simplicity we present in the follow- 
ing our general results for the case of a single decoherence 
source. The results in the case of several noise sources 
are a straightforward extension and will be discussed at 
the end. 

Finally, we can include in the model the description of 
post-selection, as required by the WV measurement pro- 
tocol. As pointed out in Ref. [THl a second QPC strongly 
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measuring the charge on the dot at any time after the 
weak measurement can effectively realize a post-selection 
in any given qubit state, I*/)- Without loss of general- 
ity we therefore consider the situation where the post- 
selection takes place immediately after the weak mea- 
surement. Within our model the post-selection into the 
state \^ f) is described by the action of the correspond- 
ing projection operator 11/ acting at the end of the weak 
measurement. 



III. GENERAL EXPRESSION FOR THE WEAK 
VALUE IN THE PRESENCE OF DECOHERENCE 

The weak value protocol we are interested in con- 
sists of preparing the double dot in a given initial state 
|\&(io)) at time to, making the system interact with 
the detector for a time r, and finally strongly measur- 
ing it in the post-selected state in |^/). The quan- 
tity of interest is the weak value of the electron's occu- 
pancy in the double dot, i.e. the value of o z {t) con- 
ditional to a positive post-selected outcome, that we 
denote by n f ( cr z( T ))*(t )- I n fact, such a quantity is 
inferred from the post-selected output of the detector, 
i.e. the average current in the QPC conditional to the 
post-selection n ( (7(r)) W ( to ) = e/r nf . (n(r))^ (to) through 
n>*(T)}* (to) =n f (T(r))^ (to) -eD/2/i)/(7 + -T_). The 
average conditional (post-selected) value of the number of 
electrons, n f (^(T))*(t ), having passed through the QPC 
during the measurement time r is 



(9) 



In Eq. ([9]) m indicates the total number of electrons 
having reached the collector and P(m \ Hf) is the condi- 
tional probability that m electrons have been transmitted 
through the QPC given that the qubit is finally found to 
be in a state represented by the projection operator II/. 
Note that Eq. Q is valid for any strength of the mea- 
surement. We will keep our analysis valid for a general 
measurement strength until specified differently. 

The conditional probabilities in Eq. can be di- 

rectly expressed in terms of the total density matrix 
R(t) of the qubit-detector system. Following the formal- 
ism of GurvitzSS and Korotkov^S a pure state of the 
qubit-detector system is described by a wave function 
= (|*t(t)),|^(t))), where a =t,| labels the 
eigenstates of ct 2 and \^ a (t)) is a many-body state of the 
QPC 



l<fj,i;r>/j, r 
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|o), 

(10) 



and |&^(i)| 2 describes the probability of finding the en- 
tire qubit-detector system in the corresponding state de- 
scribed by the creation and annihilation operators with /, 
r labeling the single particle states in the left, right reser- 
voirs, respectively. The corresponding qubit-detector 
density matrix R(t) has components 



R<rcr'(t) 



\ 



R 



R 



(m.n) 



(i) 



a matrix 
infinitely 



of 



^) r {m) 



(11) 

whose 
many 
and 



the entries 
indicates the 



Here each entry R atT ) (t) is 
dimensions are given by the 
states labelled by l,r,l' ,r' ,1" ,r" . 
l,r,l',r',l",r" ...jW.rW. Each 

^Zr 1 ...r^)t ) .ir 1 .., n (t) (the 

complex conjugate) in R^; n \t) characterizes the co- 
herences between all the possible states with m and 
n electrons detected in the collector of the QPC at 
time t. In particular, the trace of each diagonal matrix 
R^', m \t) identifies the probability that exactly m 
electrons have passed through the detector until time t, 
namely 

P(m) = tr sys [n^(t)} = tr sys {tr dot , 

(12) 

where we introduced the quantity lZ™,(t) = 
tr dct {i?^™;"^}. Also the reduced density matrix of 
the dot, p(t) = trd c t {R(t)}, can be written as 

P w = E^ (m) = E p M-P (m) w, (13) 

rn m 

where p^(t) = (t) / {tr {U^ n \t)}) describes the 

state of the qubit where m electrons have reached the 
collector. 

Besides the inherent quantum-mechanical fluctuations, 
the stochastic parameter assumes different values at 
each replica of the experiment according to its proba- 
bility distribution. In order to properly take into ac- 
count the average over fluctuations, we can first rewrite 
the conditional probability Eq. (l9| using Bayes' theorem 
as P(m\U f ) = P{m)P(Jlf\m)r£ 11} P{Il f ). The weak 
value of each run of the experiment is now weighted with 
the probability of a successful post-selection in the corre- 
sponding run of the experiment, which also depends on 
the specific noise realization. This means that the aver- 
age over the fluctuations £(t) in the weak value is properly 
taken into account by separately averaging over £ both 
the conditional average value of m and the post-selection 
probabilities. This leads to 



n f 



M r ))*(t ) = 



(Em 77 * P^f\{m\£,))P{m\Q) i 

(Em^(n / IH0)i , H0), 



(14) 



Identifying the emerging prob abilities in terms of the ma- 
trix similar to Eq. ( 12 ) yields the conditional cur- 
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rent 

(I(r)) _e E m mtr{n r (^W(r)),} 

As already pointed out the weak measurement regime 
is obtained when (5Q/Q) 2 t m < 1/-D(1 - D/Vd). The 
weak value can be identified from the coefficients in the 
expansion 



sn 



n ( = (I(r)) + nf <J(r)>^ — +0 



(16) 

We will discuss later on the validity of this expansion. 
From the definition of the weak measurement regime, we 
expect to be sensitive to coherence effects for timescales 
t m < 1/D(1 - D/V d ){il/Sn) 2 . The weak value is com- 
pletely determined by the knowledge of the probabilities 
P(m) and the conditional reduced dot's density matrix 
(p^ m H T ))£- Further analysis will now focus on the evalu- 
ation of (7^ m )(r)) 4 . 



IV. SYSTEM-DETECTOR DYNAMICS IN 
PRESENCE OF DECOHERENCE 

A. Coherent detection 

First we consider the case where a single readout by 
the external pointer takes place at the end of the mea- 
surement process, i.e. at t = r, immediately followed by 
the post-selection. In this case the coherent evolution of 
the system and detector is not disturbed by the pointer, 
and we thus name it coherent detection. 

The weak value is fully determined by the knowledge 
of the averag ed matrices (7^ m > (r)) e . We can derive 
a differential equation for (TZ^^t))^ starting from the 
von Neumann equation for the qubit-detector system, 
idtR(t) = [H(t), R(t)]. After inserting the Hamiltonian 
in Eq. (JlJ with the specific choice of a single fluctuation 
source, i.e. e(t) = e + k z ■ A(i) = A + k x ■ one 
obtains for the qubit-detector evolution 

d t R(t) = -i [ M, R(t) ] - i £(t) [ N, R(t) ] (17) 

where dt = d/dt. Here, 

M = i Aff+ + A • a x + - (Aff- + 2 e ) ® 0", , 



A = CT-k. 



(18) 



By introducing an interaction picture with respect to M 
which transforms the arbitrary operator A as Aj(t) — 



-iMt 



A 



transfers Eq. (17) into d t Ri(t) 



-i£(t) [iVj(t), Ri(t)]. Iteratively solving this equation, 



one obtains after taking the average with respect to £ (t) 



(iJ/(t)) 4 = (R(0)) £ - i / d*i<e(*i))* [JV,(ai), (i?z(0)) 



+ {-if / d Sl / d S2 <e( S i)e(s 2 )>« 

Jo Jo 
■ [JVjfo), <i?z(0)> € ]] 

+ H) 3 / d Sl f'ds, f 2 d S 3<£(si)£(s 2 )£(s 3 )>* 
Jo Jo Jo 

■[JVi(»i),[JVi(«2),[iV>(a3),<i2/(0)) f ]]] 
+ ... (19) 

Due to the 5-like time correlations in Eq. ([8]), the average 
can be performed exactly order by order. The so obtained 
integral equation for R(t) is more conveniently written in 
a differential form as 



dt (R(t)) £ = - 
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N, 



N, (R(t)). 



M, (R(t)) c 

(20) 

The density matrix is expressed as a linear combination 
of Pauli-matrices of the dot's space 



1 3 

(fl(*)>e = ,I>(*Wj 



(21) 



3=0 



where each Vj(t) is a matrix in the QPC Hilbert space 
with tr{-0o(t)} = 1. Substituting Eq. (21) into the av- 



eraged von-Neumann Eq. ([20]) yields a set of differential 
equations for Vj (t) which fully describes the averaged evo- 
lution of the qubit-detector system 



AH~,v z (t) 



AH+,v (t) 



2 7 (K - l)v x + k x k y Vy + k x k z v. 



2e v v 



AH 



2 { A ^'^l 



27 I k x k y v x + (k y - l)v y + k y k z v z + 2e v v 



2Av z 



AH + ,v y 



{AH-,v y } 



27 k x k z v x + kyk z v y + (k g - l)v z + 2Av. 



AH-,v 



(22) 



Eqs. ( 22 ) describe a set of infinitely many coupled differ- 



ential equations. They are a generalization to a density 
matrix of the results obtained in Ref. [23] for the simple 
case of a pure state. Note that such a generalization is es- 
sential in our case to properly treat fluctuations of £. One 
may further note that the fluctuations treated here dif- 
fer from Ref . [28] where the system's evolution for a given 
stochastic measurement output I(t) is studied. Consid- 
ering the matrix elements of R{t) between two sectors 



G 



(to, n) with to and n electrons having passed through 
the QPC, we trace out the QPC degrees of freedom to 
obtain a differential equation for 

{n im) {t))^ \ (vW(t).tr) , (23) 

where Vj(t) = tidet{Vj m ' m \t)}. The details of the cal- 
culation are presented in appendix [S] Here we report the 
resulting differential equation, which can be recast as a 
differential equation for wj"^ (t) , i.e. 

flivW(t) = (Go + G k + Gi) • vW(i) + G 2 • v^" 1 ^) , 

(24) 

with 











-2e 
2e -2A 
2A 



According to Eqs. ( 15|16 1 we can obtain the expres- 
sion for the weak value by solving the differential equa- 
tions ( [24] perturbatively in the regime Sfl <C ft- The 
details of the derivation are given in appendix [B] Wc 
highlight here that the perturbative solution is a valid 
approximation for r -C 1 / D(£l/5il) 2 exactly correspond- 
ing to the weak measurement regime. This finally yields 
the expression for the weak value 



n f 



J r dsv z (s) + \ J T dsn z (s) 
n • v(r) 



( 25 ) 

Here r is the duration time of the weak measurement and 
we effectively introduced a time-dependence in the post- 
selection operator H/(r) = 1/2 (n(r) • <r), n = n(r) is 
then the post-selected state at time t = r istantaneously 
after the measurement. In the notation in Eq. (25), v(s) 
is defined by v(s) = exp [(Go + Gk) s] ■ v(0) while n(s) = 
exp [(-Go + G k ) s] -n(0). 
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The above equation is obtained to lowest order in pfl 2 <C 
Vd, under the assumptions that: (i) the detector's tran- 
sition amplitude only weakly depends on the energies, 
i.e. ili r = ft; (ii) the densities of states in the QPC's 
collectors are constant, pi{Ei k ) = pi and p r (E Tk ) = p r ; 
(iii) at t — the energy levels of the detector are filled 
up to the Fermi-level so that to = 0, i.e. v^™^(0) = 



-> 



,v z ) ■ S n fi. As evident in Eq. (24), for D 
or <50 — > 0, respectively, the system evolves undisturbed, 
while for 7 — > our result reduces to that in Ref. [2H 



B. Continuous readout 

Opposed to a coherent detection, the regime of con- 
tinuous readout is characterized by the detector's state 
being frequently read out by the pointer. This limit 
is described by a sequence of readouts at times t — 
tk, k = 1...N, where the time interval between 
readouts At*, := tk+i — t k = At is the shortest 
timescale in the problem. We start by assuming At -C 
mm{l/uj(ri/Sn),l/j{n/5n),l/D(fl/sh)}. Each read- 
out corresponds to a collapse of the qubit-detector sys- 
tem in the sector of m k electrons having passed with the 
corresponding probability P(rrik]tk)- 

The evolution between two subsequent readouts is 
still excactly given by Eq. ( 24 ) with the modified ini- 



tial condition that precisely n k electrons have been 
read out at t = t\- so that v(")(ffc) = <5 n .n fc • 
(l,vi n " : ) {t k ),v i y lk: \t k ),vi nk) (t k )). In order to solve these 
modified differential equations it is useful to introduce 
a vector w(t) = (v(°'(i), vW(t), . . ,,vW(t), . . .) where n 



labels the number of transferred electrons so that Eq. ( 24 ) 
reads 



dt 



w(t) = (Mi + M 2 ) ■ w(t) 



(26) 



with M 



m,n) 



and M 



(m,n) 



(G + G k 

D 



^rl(A\ 5, 



<f - L (4) "m,n+l 



D/2) 8 m ^ 

>■> — (Gi + y 1(4)) <5 m>n + (G2 — y 1(4)) 8 m ,n+l- 

This differential equation is solved trivially in the limit 
At < l/oj(n/5n),l/~f(n/5Q),l/D(Q,/m) we are inter- 
ested in, by noting that Mi is a block-diagonalised matrix 
in Jordan form, with the solution 
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(t) = exp 



(G + G k -yl)(i-t fc ; 



r^ h+m) (t) = f m 'exp 



D 



(G + Gk- — l)(t - tk) 



(l(4) + (G!+ yl (4 ))/ 

((G 2 - f l(4))Ai + y(t - 40(1(4) + (Gi + y 1( 4 ))A«)J • v 



(27) 



(*fc) 



where / m = (y(t — ^fc}V" with m > 1. Within our 
approximation, in Eq. (27) probability conservation is 

ensured by tv{p^ stem ^> (t)} = E^o^V) = L Setting 
the reading timescales as the smallest in the problem, i.e. 
At <C 1/D, which corresponds to a continuous readout, 
the solution in Eq. (27) becomes 



v( n *)(t) = exp[(G + G k )(i-t fe )] 

■ (i (4) + Gi(t-t fc )) -v^fa) 

v^+ 1 )(t) = exp[(G + G k )(t-t i b)] 
■G 2 (i-i fc )-v^)(i fe ) 

v ("*+ m )(t) = , m> 2. (28) 

This is the limit where at most one electron pene- 
trates through the QPC between two subsequent read- 
outs.^ T he pr obabilities that either exactly one (quan- 
tum jump22E21) or zer0 electrons will have accumulated 
in the collector within a readout period time are defined 
by tr{A • v(" fe )(i fe )}o and tr{5 • v(" fc )(i fe )} , respectively, 
where the index {. . .}q denotes the zeroth component and 



A=l 



(4) 



(Go + G k + Gi) At + 0(At 2 ) 



B = G 2 At + O (At 



(29) 



0(At 2 ) indicates higher order terms for At <C 
min{l/w, I/7, 1/D}. The conditional number of trans- 



mitted electrons Eq. (15) can now be expressed as the 



sum over all permutations describing quantum jumps at 
all possible times 



n f 



("0))*(t ) = 



E mtr{ll f ■ (ft(" ll <-< mw) (r))^} 



. rriA—m 



(30) 

where j = 1...N and 7^/ mi '•■ ■ - mjv ) (7-) characterizes the 
qubit's density matrix weighted with the probability that 
exactly electrons have passed within each time inter- 



val Atk- With the definitions in Eqs. (29 I, the conditional 



number of transmitted electrons can be expressed as 



A' 



n, ( n ( r ))*(M 



n ■ E m E [A N ~ m ■ B m ] 

\ m— perm 



perm 



n- E m E [A N - m ■ B"\ cmi -v 

\ m— perm J 

(31) 



where v is evaluated at t = 0. N describes the total 
number of readouts and Eporm indicates the sum over 
all possible orders of A and B in the string of products 
A N ~ m B m . In the limit of N -> 00, At -> while keep- 
ing JV • At = t constant, the sum over all permutations 
in Eq. ( |31[) can be analytically evaluated as shown in 
appendix |C[ finally leading to 



n f 



(«(7")>*(t ) 



n ■ G 2 ■ e( G °+ G *+^+ G ^ ■ v(0) 

n . e (G + G k + Gi+G 2 )T . V (Q) 



(32) 

The weak value can be easily extracted by considering 
the expansion for small SO, /CI -C 1, so that 



n f (-^( r ))#(?o) 



1 



M T ))*(t ) - 



(33) 

The simultaneous presence of Gi + G2 and G k in Eq. ( 32 ) 
defines a new timescale 1/7 ~ l/D(Q,/8Cl) which de- 
scribes the extra source of decoherence emerging from 
the detector itself. Note also that Eqs. ( 32|33 ) are valid 
at any time. The derivation indeed relies on the per- 
turbative solution in appendix [B] which is valid in the 
limit At — > 0, hence the exact composition of subsequent 
evolution between different readouts holds at any time. 



RESULTS 



Eqs. ( 25|32|33 ) represent the main results of our paper. 
They express the weak values in the two limiting cases of 
coherent detection and continuous readout in terms of the 
system dynamics. Generally, they give rise to four dif- 
ferent timescales, which characterize (i) the system's dy- 
namics, 1/cj, (ii) the decoherence, 1/7, (iii) the detector 
dynamics, 1/D, and (iv) the backaction, 1/D ■ (Cl/Sfl) 2 . 
We realistically assume 1/D to be the shortest timescale 
in our problem. Consistently with our perturbative anal- 
ysis, we consider r <C 1/D ■ (fl/5fl) 2 . 

The weak value can be visualized via the motion of 
v(t) and n(<) on the Bloch-sphere for {v x ,v y ,v z ) 

< 1 



.1*3 IJlULll-OpilLIU L^JL V \ U X7 u y 5 

and n = (n x ,n y ,n z ) as depicted in Fig. [3] |v(i)| < 
characterizes the coherence of the qubit's state which is 
initially prepared in a pure quantum state with \v(t — 
0)|=1. 

We start with analyzing the case of coherent detec- 
tion, i.e. the dynamical evolution obtained in Eq. (25). 



In this regime the effect of decoherence essentially re- 
duces to the dynamics of the qubit alone in the presence 




11.4 " (|j 



^ 1 1 1 1 1 T I ' ' ' ' ' ' y/ti 

ti T 2 4ns Sns 12ns 0.4ns" 1 0.8ns" 1 1.2ns" 1 

(c) (d) 

FIG. 3. Evolution of v(s) and n(s) on the Bloch-sphere for (a) 
m^k and (b) m = k. (c) Weak value as a function of r for 
different 7; (d) weak value for different 7 at Ti = 0.10 ns and 
T2 = 1.31ns, marked in panel (c). In all plots the parameters 
are chosen as eo = 20 fieV, Ao = 3 ^ieV, k ~ — 0.89^ + 
0.45 e z , v(0) = (0.5,-0.33,0.80) and n = (-0.20,0.75,0.63). 



n/ «r, (t)>^ n/ <oi W>;™' 




4ns 8ns 12ns 



(c) 

FIG. 4. Weak value as a function of r for (a) different k (8 = 
corresponds to k | m) and (b) different n (<f> = corresponds 
to n || m); (a) decomposition into two fluctuation sources 
where (i) 7 = l.O^eV, k = + ^e z , (ii) ki = ^ej, 

k2 = ^e 2 , 71 = 775 A^V 7 = 775 M e V, (hi) ki = m, k 2 = 
mx, 71 — (k ■ m) /ieV, 7 2 = (k • mx) /xeV; all the other 
parameters are as in Fig. [3] 



of decoherence. Accordingly, the weak value presents dif- 
ferent behaviors in different regimes defined by the dura- 
tions of the measurement, r, compared to the remaining 
timescales 1/ui, 1/7. 

For long measurement durations r 3> I/7, the qubit's 
state generally relaxes towards a fully statistical mixture, 
v(i) = 0, as shown in Fig. [3^c) , except for the special 
case k \\ m, which will be discussed later. Consequently, 
at timescales where decoherence effects come to play we 
obtain max| nf (^(t))^?^| < eD and, hence, the peculiar 
characteristics of weak values are washed out. 

For measurements shorter than the decoherence time, 
the results depend on the system's evolution timescale. 
For short measurements, t <C 1/w, which correspond to 
negligible dynamics and fluctuations, the vectors v(s) 
and n(s) in Eq. (251 are constant so that a measure- 



when |1 + n • v(t)| ^> 1 which leads to r = n ■ — ± At 



with At < 



ment of the averaged detector's response trivially reflects 
the weak value of the observable a z independent from 
the measurement duration time t, nf (<j z (t — > 0))^^ = 

Re«f|*,|*)/<f|*>). 

If the decoherence is weak enough, t <C 1/7, the sys- 
tem's dynamics for intermediate durations of the mea- 
surement, t 3> 1/w, however are appreciable, then both 
vectors n(s) and v(s) precess about the eigenvector m 
of Go with a frequency of 2lj. Note that, due to its 
backward-in-time evolution, n(s) precesses in the oppo- 
site direction as compared to v(s). Due to the oscillatory 
dependence of the denominator in Eq. (251, peculiar weak 



values may occur for properly fine-tuned measurement 
duration times. As an illustrative example weak val- 
ues much larger than 1 are realized for orthogonal states 



I sin 2 a I 2 1 



where a denotes the altitude angle 



between v(0) and m. The effect of the strength of the 
fluctuations on short and intermediate measurements is 
reported in Figs. [3^c) and[3^d), resulting in a trivial de- 
cay towards the steady state value. 

As already anticipated, an important role is played by 
the "direction" of the noise term. Its effect is illustrated 
in Fig. 4ja). As long as k becomes more and more paral- 
lel to the m, the relaxation time toward the steady state 
becomes longer. In the limiting case, it relaxes to a par- 
tially coherent state. The direction of the post-selection 
orientation plays a similar role, as shown in Fig.Qb). We 
also note that, in the presence of multiple decoherence 
sources, even for given effective direction and strength 
of the fluctuating term, the decay of the weak value still 
strongly depends on the relative directions between dif- 
ferent sources (c.f. Fig. E[c)). 

The case of a continuous readout presents signifi- 
cant differences which are visualized in Fig. [5] While 
in both cases a sufficiently strong decoherence destroys 
the peculiarity of weak values, in the continuous read- 
out regime the detector itself introduces decoherence on 
a timescale of its backaction. Interestingly, the weak 
value approaches nf W z (t))^ ) -> n z + 0{5n 2 /fl 2 ) for 

l/D < 1/uj < I/7 ~ l/D(n/Sn) 2 < t independent 
from the direction of the fluctuations which corresponds 
to a measurement of the z-component of the post-selected 
qubit's state. Such long-lasting measurements, however, 
are beyond the interesting weak measurement regime. 
The dependence at timescales shorter than the deco- 
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n/«n(r)>J; 
1.2 
0.8 
0.4 



n/<CTi(T)>?eij 





- y=0.O/ieV 




y=0.5/ieV 




y=l,0/teV 



coherent 
continuous 
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(a) 



4ns 8ns 12ns 

(b) 



n/ «7- E (T)>; 



nj <o--(r>>^ i k ) 




(c) 



(d) 



FIG. 5. (a) Continuous detection for different 7 and (b) com- 
pared to the coherent readout for 7 = 1.0 fieV; D = 0.270 eV, 
SQ/Q = 0.001, other parameters as in Fig. [3] (c,d) Continu- 
ous detection and coherent readout as in (bjfor 7 = 0.1 /ieV, 
v(0) = (0.5,-0.33,0.80) and n = (-0.20,0.75,0.63), where 
D = 0.270eV in (c) and D = 0.027 eV in (d). 



herence time behaves more remarkably, where we are in- 
deed in the regime of a weak measurement. Here, the 
outcome is also strongly effected by the detector. In 
Figs. [5](b-d) the results of continuous readout and co- 
herent detection are juxtaposed. We clearly see that a 
continuous measurement regime might enhance the weak 
value as compared to a coherent readout with the same 
decoherence strength. In fact the detector not only in- 
troduces a new decoherence time, but it also affects the 
system's dynamics at r <C \/D(Vl/8Q) 2 by determining 
its precession axes. This implies that the enhancement 
effect is non-universal and depends for instance on the 
post-selection direction or the noise strength, as shown 
inFig.gb)-(d)._ 

The perturbative solution of the system's dynamics un- 
derlying the result of the coherent detection allows us to 
discuss, to some extend, the validity of the weak value 
expression in Eq. |25| As a first check we can require that 
the second order is irrelevant compared to the first or- 
der contribution discussed so far, i.e. | nf (I ( T ))^(tg) \ ' 
Stf/n 2 < I nf (I(t))^| ■ Sn/Q, leading to the condi- 
tion 



77(7-) 



sn 
IT 



Dt 
~2~ 



V x {t) n x + v v (t) n y 



\ Jo ds v z(s) + \ Jo ds %( s ) 



< 1 



(34) 

This imposes a restriction on the validity of the weak 
value's result also within the regime of weak measure- 
ment, as shown in Fig. [6] Indeed, for specific qubit's 
parameters eo and Ao and particular boundary condi- 
tions v z (0) and n z , the numerator in Eq. ( |34[ ) vanishes 
at finite duration times for r <C 1/7 so that the per- 
turbation is valid at discrete times which depend on the 



lilllllllllii.. 



4ns 8ns 12ns 

(a) 




FIG. 6. Validity perturbation; (a) left-hand-side of Eq. (341 
and (b) the corresponding weak value for 7 = 1.0/j.eV, D = 
0.270 eV, SQ/Q = 0.001; all the other parameters are as in 
Fig-i 



chosen parameters (c.f. Fig. K3J). For r » 1/7, on the 
contrary, 77 oc r 2 e~ CT where c > and the perturbation 
is asymptotically valid unless k = m = e^.. 



VI. CONCLUSIONS 

In this work we have addressed the effects of deco- 
herence on weak value measurements involving post- 
selection. We have considered the paradigm model of 
a charge measurement in a double quantum dot by a 
nearby quantum point contact where we have included 
fluctuations of the parameters due to external noise 
sources. After deriving a general expression for the post- 
selected signal (current) in the QPC in terms of the re- 
duced density matrix of the qubit, we have evaluated it 
explicitly in two different regimes set by the detector's 
readout, namely continuous vs. single-time detector's 
readout. 

We have characterized the weak value's dependence on 
the various parameters of the system. In particular, we 
have shown that statistical fluctuation's of the qubit's pa- 
rameters generally reduce the weak value into the clas- 
sical range for measurements longer than the decoher- 
ence time. On shorter timescales we have determined 
a boundary for the region of validity of the weak value 
result. Remarkably, there, due to its interplay with the 
system-dynamics, the continuous readout can lead to an 
enhancement of peculiar weak values as compared to the 
case of coherent detection. 
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Appendix A DERIVATION OF AVERAGED 
RATE EQUATIONS 

In this appendix we derive the differential equa- 
tion (24) out of Eqs. (22) of the main text. In the fol- 



lowing, the derivation is presented exemplarily only for 
v x — tr{?)i m ' m ' ) (t)} since the other terms {v y m \ ) 
are treated completely analogous. It is useful to perfom 
a "Laplace" transform for the whole matrices 



s (m,n) 



lim 

5^0 



- (m,n) 



(t) exp [i (E + i 5) t] dt (A.l) 



in order to include the initial conditions of the differential 
equations that to = electrons have penetrated through 
the collector at t = and, hence, the qubit-detector sys- 
tem is in a pure state. Here, 5 > ensures the conver- 



gence of the integral. A high energy cut-off is introduced 
concerning the inverse transform, so that the upper limit 
of the integral in the inverse transform is A — > oo 

We can write the differential equations for the Laplace 
transformed components, v^ m ' m \t). To this regard, the 



matrix products in Eq. ( |22| ) can be easily calculated by 
evaluating each (ro,n)-block as introduced in Eq. (11). 
Due to the fact that the definition of AH + includes terms 
proportional to a], a r and aj a\ , diagonal blocks of AH + 
are non-zero whereas the diagonal blocks AH~ ( m < m ) 
vanish. Moreover, since AH^ 1 consists of combinations 
a\ a r +h.c. which raise or lower, respectively, the number 
of electrons in the detector by exactly one electron, only 
the off-diagonal blocks AH ± an d AH ± ( m +^ m ) 

neighboring the diagonal blocks are non-zero. Conse- 
quently, a block-wise evaluation of the matrix products 
can be written as 



■ Air) ' =tf m ' n ~ 



( A -2) 

(vj ■ Aff+) {m,V) = v( m ' n -V . AH+ + vf l ' n) ■ AH+ ^ + v^ m ' n+1) ■ AH+ {n+1 ^ (A.3) 



and analogous for (AH + / ■ Vj)( m ' n \ Note that v{ m ' k ^ ■ AH ( fc ' ra ) still describes a product of matrices which are of 



infinite dimension. Thus, applying the "Laplace" transform to Eqs. (22) for j = x by considering exemplarily k = e 2 
leads to 



tr 



= -£- 



tr 



£~ l 



+ 



tr i -2 7 vr ,m) (E) - 2 e (E) 



=«a(-E) 



=«i(£) 



£^jj+(m,m-l) _ ~(m-l,m)/m _|_ ^jj+(m,m) _ ~(m,m) /g\ _|_ ^Jj+ (m,m+l) _ -(m+l,m) 



=«c(S) 



=v d (E) 



fi(m,m-l) , ^jj+(m-l,m) _|_ £j(ro,m) /£A . (m,m) _|_ ~(m,m+l) /£\ , ^jy+(m+l,m) 



AJJ- (m,m-l) , ~(ro-l,m) j-^-j _|_ ^jf - (m,m+l) . -(m+l,m) 



V 



=%(-E) 



+ W^ m ' m_1 ^(i?) • AH~ ( m_1 ' m ) + {j( m ' m+1 ) (£?) • Aii" - (m+l,m) 



=S fc (E) 



(A.4) 



where ,£ 1 denotes the inverse " Laplace" transform. We analyze the various terms separately. We start with observing 
that £~ l [v a (E)] = — 2jtr{v^ n ' ,m \t)} and £~ 1 [vb(E)] — — 2 e tr{vy m ' m ^ (t)}. Moreover, due to the cyclic invariance 
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of the trace £ 1 [id(E)] + £ 1 [v g (E)] — 0. Employing the explicit expression of AH ± , reduces the evaluation of the 
remaining terms to the calculation of 



(Ajy+Cm.m+i) .-(m+i,n)\ = y (-l) m 2r» / +ir + ■ (i 



^ (m+l,m) \ 

■lmr m );{l[r[...l'T' n ) , " ' ■"■-r-'-.-r- V- 7 / (/ in .. .i m r m i m+ ir m+ i) ;(iiri . . .l' n r', 



(A.5) 



and 

( A jj+ (m,m-l) . fi (m-l,n)\ = _ 2 fi , / fi (m-l,m)\ 

V 3 /(/iri...i m rm)i(2iri-inO ^ ' 

■ r> n /" «(m— 1,7 



(hr2---l m r m );(l' 1 r' 1 ...l' n r' n ) 
(hr 1 l 3 r 3 ...l m r m );(l' 1 r' 1 ...l' rl r' n ) 
(Iir 1 l2r 2 l4r4---l m r m )-(l' 1 r' 1 ...l' n r' n ) 



9 /V,( m_1 > m ) 

+ (-i) m 2n lmrm ■ U m ~ l > m A (a.6) 

V 3 Hhri-l m -ir m -i);(lir' 1 ...l' n r' n ) 



for products concerning the off-diagonal boxes and 



( AH+ -v[ m ' n) ) = 2 (E ri +...+ E r -E l -...-E l )■( v ( ^ n) \ 

V 3 J (l 1 r 1 ...l m r m );(l' 1 r' 1 ...l' r' ) ml / (hri ...Z m r m );(Ziri ...Z' r' ) 

(A.7) 

for products with non-zero matrices on the diagonal (m, m). Here, the indices {l\r\ . . . l m r m ); (^r^ . . . l' n r' n ) precisely 
determine the scalar entries of the matrices. The analogous multiplications involving AH~ are evaluated by replacing 
Q — > SQ. In principle, these products are not restricted to the usage of Vj(t) and, thus, are valid for a generic matrix 
of the dimension of £}j m±1,m ), it is sufficient to consider products where AH^ 1 is on the left, since in the end we need 

to evaluate the traces and all products can always be arranged such that AH^ is on the left by relying on the cyclic 
invariance of the traces. 

To lowest order in pfl 2 <C Vg, we find that 



(^(m,n)\ f ~(rn,n)\ 

3 J (hr 1 ...l k r k ...l m r m );(l' 1 r' 1 ...l' n r' n ) 2 ^ 3 hhri...l 



^2 V _ Hhri...l k r k ...l m r m );(l' 1 r' 1 ...l'„r'„) ^ ^ 2 V _ ' {hri...l m r m );{l' 1 r' 1 ...l' k r' k ...l' n r' n ) ^ . . 

fa E + E lh -E r „+i (8 + 2<y) ~ ' fa E + E Vh - + i (S + 2 7) " " 

if the sum runs over indices occuring in v^ m ' m \E) with the abbreviation E = E — E^ — ... — Ei k l — Ei h+1 — . . . — 
Ei m + E ri + . . . + E rk l + E Tk+1 + . . . + E Tm . This result holds within the approximation where the energy-levels 
of the detector's reservoirs are almost continuous so that the sum can be replaced by an integral, i.e. J2i k r k 
/ pi(Ei k ) p r (E rk ) dEi k dE rk . Additionally, the hopping amplitudes are assumed to depend weakly on the energy 
levels, hence, being constant, i.e. f2; r = f2, and also the density matrices of the emitter and collector, respectivley, are 
approximated to be constant, with pi(Ei k ) = pi and p r (E rk ) = p r . In order to understand that the integral vanishes, 
it is helpful to realize that the all entries of t)j m ' n - ) essentially characterize higher order retarded Green's functions 
which describe the averaged evolution of the density matrix. These Green's functions have poles in the lower half 
of the complex plane proportional to E + Ei k — E Tk + iS]^ 1 which can be shown by iteratively solving the Laplace 
transformed differential equations (22) for each v^ n ' n \E). Thus, a contour integral yields zero since the integrand 



decreases oc 1/E 2 . On the contrary, if the sum does not include a summation over indices of Vj' n \E), it can be 
calculated within the same approximation as being 

17 2 V^ 1 -^-iir Pl p r n 2 (V d + E) n- l -D (A.9) 

to lowest order in p J7 2 <C Vd with D = 2 it pi p r SI 2 Vj^- Concerning this, the integral is separated into a singular 
part and the principal value part. While the principal part redefines the energy-levels, the singular parts leads to the 
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presented result by relying on the Sokhatsky-Weierstrass theorem. Thus, one obtains 
V C (E) « j B • ( 2 tr i) (E)} + 2 i f tr {v(m-i,m-i) {E) ^j 

-tr|4 m_1,m) (0)- (iiE-6-2j)H m ' m) + % - Aff +(m ' m) J • Afl" +(m - m - 1) | , (A.10) 
v e (E)K-~D. (2 tr (£) } - 2 * ^ tr {« (£) }) 

- tr |a//+(™.™+1) . ^(t^- 5- 27) i^ + Mn+l) + * AJ y+(m+l,m+l)^ . fi (m+l,m) (q) | > ( AU ) 

W « 4 D ■ (- 2 ^ tr {v^-^(E)} + tr ^(^O}) 

- tr|^ m-1 ' m >(0) ■ ^(i £7 — <5 — 2 -y) l< m ' m ' + % - AH +( - m ' m ^ ■ AH~ ( m < m ~V 1 (A.12) 
« ^ • (2 f tr + 2* ^ tr {e' m) (^)}) 

-tr|Aff-( m ' m+1 > • ^(i E 1 - <5 - 2 7) i(" l + 1 . m + 1 ) _ i A# + ( m + 1 ." l + 1 )^ . 7 }( m + 1 ^)(0)| . (A.13) 

which couples (m,m)-blocks to (to — l,m — 1)- and (to + l,m + l)-blocks, respectively. Note that it is trivial to take the 
inverses since the corresponding matrices are diagonal. Similar calculations eventuate in analogous expression for all 
the other terms by replacing m by to + 1 or m — 1, respectively, i.e. if(E) = v e (E)(m — > m— 1), vu{E) = v c (E){m — > 
to + 1), Vk{E ) = Vj(E)(m — > m — 1), #/(.E) = Vi(E)(m — >• to + 1). Inserting the achieved results for v a {E) to 

into Eq. (A.4| and iteratively solving in to, yields the desired differential equation for i^"^(t) := tr Iwj" 1 '" 1 ^)! with 

(A.14) 



j = x after performing the re-Lapclace transform, i.e 

12 



*£">(*) = (-27- f (l + ^)) + f (l-'gf) ^- 1) (*)-2e ^(t). 

Likewise calculations for each vj (i) finally end up in the rate equations 

) vt\t) + f (l + ^) ,r X) (*) - C f ^ (*) + U f (A.15) 





(*) = 






.(m 
w 




+ w 




(*) = 








(-27- 






(*) = 








-?0 


+ 7^ 



(A.16) 



*4 m) (*) + ~ 



(A.17) 

for the special case of considering fluctuations in the z-direction. The general case is treated analogously. 

Appendix B PERTURB ATIVE SOLUTION OF THE RATE EQUATIONS 



In this appendix we present the perturbative solution for 5fl <C Q of the rate equations ( 24 ) which describe the 



averaged evolution of the density matrix with the initial conditions that to — electrons have penetrated to the right 
reservoir at t — so that the initial state of the system is described by Vj m \t = 0) = S m fi ■ v< f\t = 0), j = 0,x,y, z. 
Hereto, we introduce the Fourier transform 



13 



with respect to m. Hence, (TZ 1 - 171 ^^))^ becomes (TZ(q,t))^ — 1/2 • (v(q,t) ■ a) which eventuates in an expression of the 
conditional current in terms of Vj(q,t) by comparing to Eq. (15), i.e. 



n ( ( / ( T ))*(to) = ~ — d 1 ln 



tr \ ll f ■ (v(q,t) ■ a) 



(B.2) 



q=0 



If assuming that v(q,i) is an analytic function of q and t, it can perturbatively be expanded in order to describe the 
evolution of the system's density matrix within a weak measurement regime 



V M) = 2^ Un(?,*) I 7^ ) 



(B.3) 



where n denotes the order of perturbation. Substituting this perturbation into Eq. (B.2 1 the zeroth order contribution 
reads 



(Hr)) 



■i— d„ ln 

T 



tr \n f ■ (u {q,t) ■ a) 



(B.4) 



q=0 



and the first order contribution, i.e. the averaged weak value for the current, is identified as 

weak _ • 



n f U( r ))*(t 



tr |n / .(u o (? = > t).ff)|-tr ^ILfid^iq^l^-a 
tr {n / -ui(g = 0,t)-<7)j-tr \u f ■ (d q u (q,t)\ q=0 ■ a) 



tr \ n f ■ (u (q = 0,t) ■ a) 



(B.5) 



Thus, the weak value is completely expressed in terms of the averaged density matrix. Further analysis will now focus 
on the evaluation of u n (q, t) for n = 0, 1, aiming at finding an illustrative expression for the weak value. Inserting the 
power series of Eq. (B.3 1 into the Fourier transformed rate equations (24) leads to a set of differential equation for 
each u n (q, t). The resulting equations for the lowest orders are 



d_ 
dt 

d_ 
dt 

d_ 
dt 



Here, we have introduced 



u (g,i) = 


Go - 


FG k - 




1(4) 


• u (g,t), 




ui(g,t) = 


Go - 


\-Gk — 


f (l-e'«) 


1(4) 


• ui(g,t) - 


D(l-e^)-G i-u (q,t) 


u 2 (q,t) = 


Go - 




f (l-e") 


1(4) 


• u 2 (g,i) - 


~2 'G g - u (q,t) 



-L>(l-e i9 )-G i-exp 



fGo + Gk-y(l-e <9 )-l w )< 



Ui (£?,*) 



(B.6) 
(B.7) 

(B.8) 



Gqi 



1 





10 



G 



10 



/ 
10 
10 

Vo o o o 



and G„ 



/l-e iq \ 

l + e*« 

l + e l « 

\ l-e iq j 



(B.9) 



Higher orders in u„(g,i) are not relevant for the expression for the weak value. Pertinent for the expression of the 
weak value are solutions for the special cases u n (q = 0,t) and d q u n (q, f)| q =o- 

The initial conditions are equivalent to u a (q,t = 0) = (l,Vx(t — 0),Vy\t — 0),vi (t — 0)). Thus, u (q,< = 0) 
does not depend on q initially and u n (q,t = 0) = for n > 1. Furthermore, d q u n (q,t — 0)| 9= o = for all n € N and 
tr{(p(t))j} = 1 implies (u (q — 0,t)) = 1 and (u n (q = 0, t))o = at any time t. Additionally, d t u n (q = 0,t) = at 
any t to keep tr{(/9(t))g} unchanged. 

The perturbative solution is obtained iteratively with v(t) = exp [(Go + Gk) i] ■ v(0) where v(0) = 
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(«i 0) (o),4 0) (Q).«* 0) (o))- lt reads 

u o(ff = 0,t)=( v ( t) ), 11^ = 0,*) «^ = 0^)=( i)tMt )e x + % (i)e,) 



1 

v(t) 



9 g ui(g,i)| g= o = iD 



d g u 2 (q,i)| 9=0 = 

([G + Gk]" 1 ■ exp [(Go + G k ) t] - [G + Gk]" 1 ) ■ v(i - 0) 
([Go + Gk]" 1 • exp [(Go + G k ) t] - [G + G k ] _1 ) • e 2 

Herefrom we conclude that 

tr{n, • (uo(? = 0,t) -<j)} = l + v(r) • n, tr {il/ • (m(g = 0,t) • <r)} = 

and 



(B.10) 
(B.ll) 

(B-12) 
(B.13) 



t T 

:tr{n / -(S 9 u 1 (g 1 t)| g=0 .<r)} = -ty ^tr{n / -(S,ui(g,s)| 9=0 .a-)}ds = D J v z {s) + n z (s) ds . (B.14) 



Inserting Eq. (B.14| and Eq. (B.13) into the expression for the conditional value in Eq. (B.2) then yields the expression 
for the weak value, i.e. Eq. (25), which is presented in the main text. 



Similar, the second order contribution to the weak value is evaluated by noting that 



itrjlL/ ■ (d q u {q,t)\ q=0 ■ tr) I 



Dt 



1 + v(t) • n 



and 



tr{n / .(u 2 (g = J *)-(T)} 
which eventuates in the expression 



Dt 



v x {r) n x + v y (r) n y 



,j, ^wcak(2) _ _ eD 2 T V x (t) n x +V y (T) 7l y 

ntVW'vfa) ~ 2 ' 1 + n-vM ■ 



(B.15) 



(B.16) 



(B.17) 



Note that the second order contribution in Eq. (B.17) has the same characteristics, i.e. the same denominator, as the 
first order term in Eq. (25) which implies analogous conditions for divergencies or peculiar weak values. 



Appendix C THE SUM OVER ALL PERMUTATIONS 



Here, we show how the sum over permutations in Eq. (32) is performed. Defining the auxiliary function 
f(m,At,N) = ^ perm [A N - 



1 B m \ , the numerator is evaluated as follows 

J perm J 



N 



n - I ^ m - f(m,At,N) • v 



\m=0 



in ■ d„ 



in ■ d„ 



N 



-imq 



^2 /(to, At, N) i 

N 

E E [A N - m -(Be-^r] pe 



jn— perm 



in • d q [(1 + (G + G k + Gi + G 2 e - <9 )At) 



N' 
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where the last step is valid since J2m=o Sperm [A N m • B m ] grm = (A + B) N and due to the definitions in Eqs. ( 29 1 
In the limit N ->• oo, At -t 0, N ■ At = r, this readily yields rn ■ G 2 ■ exp((G + G k + Gi + G 2 )r) • v(0) 
denominator is treated analogously. 
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